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Abstract 

The flow of fluid confined between a heated rotating cylinder and a cooled 
stationary cylinder is a canonical experiment for the study of heat transfer 
in engineering. The theoretical treatment of this system is greatly simplified 
if the cylinders are assumed to be of infinite length or periodic in the axial 
direction. In these cases heat transfer in the laminar regime occurs only through 
conduction as in a solid. We here investigate numerically heat transfer and the 
onset of turbulence in such flows by using both periodic and no-slip boundary 
conditions in the axial direction. The influence of the geometric parameters is 
comprehensively studied by varying the radius ratio (0.1 < 77 < 0.99) and the 
length-to-gap aspect ratio (5 < F < 80). Similarly, a wide range of Prandtl, 
Rayleigh, and Reynolds numbers is explored (0.01 < a < 100, Ra < 30000, 
and Re < 1000, respectively). We obtain a simple criterion, Ra < a(? 7 )F, 
which determines whether the infinite-cylinder assumption can be employed. 
The coefficient a is well approximated by a cubic fit over the whole ry-range. 
Noteworthy the criterion is independent of the Prandtl number and appears 
robust with respect to Reynolds number even beyond the laminar regime. 
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1. Introduction 


Instabilities driven by the combination of rotation and thermal gradients 
determine the dynamics of complex geophysical, astrophysical and industrial 
flows. Simple models of such flows can be tested in laboratory experiments of a 
laterally heated differentially rotating annulus, for which extensive information 
about the physical mechanisms and flow regimes can be found in the literature 
(see [1] for a review on this topic). The case of rotating heated inner cylinder 
and stationary cooled outer cylinder (RHISCO) is a model for the cooling of 
rotating machinery, the solidification of pure metals, techniques of chemical 
vapor deposition, rotating-tube heat exchangers and nuclear reactor fuel rods [2 
ElllllS]. The geometry of such an experimental apparatus is fully specified by 
the length-to-gap aspect ratio T = h /{tq — Ti), and the radius-ratio rj = rijro^ 
where ri and Tq are the radii of the inner and outer cylinders, and h is their 
height. 

Depending on the geometry, RHISCO experiments used in the literature can 
be classified in two groups. The first group of experiments [SlEli is charac¬ 
terized by long cylinders T = h/{ro — Vi) > 100 and narrow gap rj = Tijro < 1 . 
Ali and Weidman [9] performed a detailed linear stability analysis of such flows 
using axial periodicity and reported on the influence of the Prandtl number 
(cr) and rj on the stability boundaries. Their results showed a good agreement 
with [ 6 ] and, to a lesser extent with [?]• Ali and Weidman attributed the dis¬ 
crepancies to the limitations of linear stability theory and the infinite-cylinder 
idealisation to capture the experimental details. A similar linear stability anal¬ 
ysis dU] reported good agreement between numerical and related experimental 
results [ 8 ]. Nonlinear simulations for small temperature gradients were provided 
by Kedia et al. [la who quantified the heat transfer across the system. A second 
group of experiments embraces apparatuses with moderate aspect ratio and wide 
gap. Ball and Farouk [EIODj] reported heat transfer measurements as well as 
the sequence of flow transitions using an experimental setup with T = 31.5 and 
77 ^ 0.5. Subsequent numerical simulations m for r = 10 and 77 = 0.5 provided 
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insight on the bifurcations structure of the system. However, the results showed 
significant discrepancies with experiments suggesting strong effects of the axial 
boundaries. 

An accurate numerical simulation of the axial (Ekman) boundary layers in 
flows between long cylinders entails a substantial computational cost, especially 
for rapid rotation and large temperature gradients. The assumption of axial 
periodicity reduces the computational effort because the Ekman layers are not 
present and only a short central fraction of the apparatus needs to be simulated 
to estimate transport properties. Moreover, under this assumption variables can 
be expanded as a Fourier series in the axial direction, which greatly simplifies 
the numerical approach and enables the use of more efficient solvers. 

In this paper we determine under what conditions periodic boundary con¬ 
ditions can be employed to describe the dynamics of RHISCO experiments. In 
particular, we provide criteria to distinguish flow features that arise from the 
interplay between differential rotation and temperature, from those which are 
mainly determined by the axial boundaries or end walls. We compare the flow 
dynamics by using both physical (no-slip) and axially periodic boundary con¬ 
ditions in our numerical simulations. We show that axial periodicity renders a 
good approximation of laboratory flows as long as the Rayleigh number Ra is 
small. We provide a simple criterion that determines whether heat transfer in 
the laminar flow is conductive or convective. In particular, conductive profiles, 
which enable the use of axially periodic boundary conditions, are obtained as 
long as Ra < a{r])T + b{rj). In addition, the coefficient b may be neglected for 
r > 15. 

2. Specification of the system and numerical methods 

We consider the motion of an incompressible fluid of kinematic viscosity v 
and thermal diffusivity k confined in the annular gap between two rigid and 
concentric rotating cylinders of radii and Tq. The inner cylinder rotates 
at constant angular velocity whereas the outer cylinder is kept at rest. A 
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Figure 1: (color online) Sketches of geometry and boundary conditions for the 
two cases considered in this paper, (a) Stationary insulating endwalls (no-slip 
boundary condition for the velocity and zero flux for the temperature) and (b) 
axially periodic boundary conditions. The temperature profile is superimposed 
on the right hand side of each figure. 

radial thermal gradient is considered by setting the inner and outer cylinder 
temperature to = T + ATj2 and To = T — AT/2, respectively. Here T 
is the mean temperature of the fluid. The axis of the cylinders is vertical, 
i. e. parallel to the gravitational acceleration g. We study flows with stationary 
end walls and with axially periodic boundary conditions (Figures (a) and (6) 
respectively). The latter model the case of infinitely long cylinders, whereas the 
former reproduce experimental boundary conditions (no-slip for the velocity and 
thermally insulating end walls for the temperature). 

2.1. Governing equations 

We consider the Boussinesq approximation including centrifugal buoyancy 
in an inertial reference frame as described in |16j . The dimensionless governing 
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Name 

Symbol 

Formula 

Operation range 

Grashof number 

G 

agNTd?/v^ 

0 < G < 12000 

Relative density variation 

e 

aAT = Apjpo 

0 < e < 0.025 

Prandtl number 

a 

V ! K 

0.01 < cr < 100 

Aspect ratio 

r 

hjd 

5 < F < 80 

Radius ratio 

V 

Til To 

0.1 < ?7 < 0.99 

Reynolds number 

Re 

Q.rid/u 

0 < Re < 1000 


Table 1: Dimensionless parameters. Here k and a are respectively the thermal 
diffusivity and the coefficient of volume expansion of the fluid, Ap is the density 
variation associated with a temperature change of AT. Last column indicates 
the range of values covered by the simulations shown in this paper. 

equations are 

(i9t + V • V)v = —Vp + V^v + GTz + eTv • Vv, (la) 

((9t + V • V)T = ct-^V^T, (lb) 

V-v = 0, (Ic) 

where v = (m, u, w) denotes the velocity field vector and T is the deviation of 
the temperature with respect to T. The length, time, temperature and pressure 
scales chosen to make the set of equations dimensionless are the gap width 
d = To — Ti, the viscous time jv, the temperature difference between the 
cylinders AT and (vjd)'^ respectively. There are six independent dimensionless 
numbers (see table [^. The term eTv • Vv accounts for centrifugal buoyancy, 
including secondary effects stemming from differential rotation or strong internal 
vorticity [16j . The equations are solved in cylindrical coordinates {r,9,z) 

2.2. Numerical methods 

In the axially periodic case the onset of instabilities was determined via linear 
stability analysis of the basic flow as in [16) . Fully nonlinear simulations were 
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performed using the Boussinesq-approximation m, which was added with the 
heat equation to the hnite-difference-Fourier-Galerkin (hybrid MPI-OpenMP) 
code of Shi et al. m- A time-step (5t = 2 x 10 ® viscous time units was used in 
all computations. 

For rigid end walls the governing equations were solved using a second-order 
time-splitting method. A pseudo-spectral formulation is used for the spatial 
discretisation, with the Fourier-Galerkin method in the azimuthal coordinate 9 
and Ghebyshev collocation in r and z. The code is based on a previous hydrody¬ 
namic code |18) . which has been extended with the Boussinesq-approximation 
of [TB] and parallelised as in m- Details can be found in J. M. Lopez’s PhD 
thesis m- The numerical resolution has been chosen to ensure that the inh- 
nite norm of the spectral coefficients decays at least four orders in magnitude. 
Time steps as small as = 1 x 10“^ viscous time units have been required for 
numerical stability and accuracy of the second-order temporal scheme. 


3. Conductive and convective basic flows 


The assumption of axial periodicity allows to considerably simplify the cal¬ 
culation of the basic flow. The radial velocity is zero and the rest of variables 
only depend on the radial component. Under these conditions an analytical 
solution can be found by imposing a zero axial mass flux. This reads 


B 


Vb{r) = Ar-\ -, 

r 

Wb{r) = G(cir^ - r^) + (c{rl - r^) + 
ln(r/ro) 1 


Tb{r) = 


Inrj 



(2a) 

(2b) 

(2c) 


The azimuthal velocity Vb is the classical Gouette flow, Tb corresponds to the 
temperature in a conductive regime and Wb is the axial velocity profile induced 
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Figure 2: Radial profiles of (a) axial velocity (computed for G = 2000) and ( 6 ) 
temperature as a function of the radius ratio 77 . 


by gravitational buoyancy. The parameters A, B and C are 


Rco — rjRe Re — rjReo 

1 + 77 ’ ^( 1 - 77 )( 1 -? 72 )’ 

41n 77 + (1 — 77^)(3 — 77^) 

16(1 — 77^)((1 + 772)ln77 + 1 — 77 ^) ’ 


(3) 

(4) 


and the non-dimensional radii are = 77/(1 — 77 ), Tq = 1/(1 — 77 ). 


Note that the conductive temperature (2c) only depends on the radial ge¬ 
ometry. Figure [^illustrates the variation of the axial velocity and temperature 
with 77 . 


3.1. End walls effect: influence of fluid properties and geometry. 

The presence of end walls in the system leads to non-zero radial velocities 
close to the end walls, commonly referred to as Ekman layers, and enables 
convective heat transfer across the gap even for the basic flow [201 mi- This 
switches on the contribution of terms of the heat equation, which were zero 
under the axial periodicity assumption. The temperature equation reads 

dtTb + UbdrTb + -deTb + Wbd.n = (5a) 

r 
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In the infinite case the steady basic fiow has no radial velocity {ub = 0) and 
the temperature (Tf,) depends only on r, so all terms in the left hand side 
are zero. Consequently, the basic flow is unaffected by changes in the fluid 
properties. In contrast, the influence of a may be expected to play a significant 
role in bounded systems. It modifies the basic flow with respect to the idealised 
periodic situation and may become an important source of discrepancies between 
both systems. 

Figurej^a) shows colormaps of the temperature in a meridional section (r, z) 
for rj = 0.5, Re = 30, G = 2000, F = 10, and a = 0.1,1,10. It clearly illustrates 
how changes in a alter the temperature of the basic fiow in a finite system. 
Radial and axial temperature profiles are shown respectively in figures (b) at 
mid-height and (c) at mid-gap as a function of a. For ct < 1 (gases or liquid 
metals), and sufficiently far from the end walls, the temperature profiles are 
nearly independent of z and thus deviations from the conductive regime (cr = 0) 
are negligible. The axial dependence becomes stronger with increasing a, re¬ 
sulting in radial temperature profiles that differ substantially from a conductive 
profile even at mid-height. This temperature regime is usually referred to as 
convective [52J US] , and is characterised by sharp thermal boundary layers and 
enhanced heat transfer. 

The heat transfer across the gap is here characterised by the Nusselt of the 
inner cylinder 

I-r%drT\r,dzde 

< Nu, >t= ^- ( 6 ) 

Qcond 

where < . >t indicates time average and gqcond denotes the dimensionless heat 
flux of the laminar (conductive) flow in the periodic case. To simplify the 
notation the subindex i is omitted henceforth. Note that the Nusselt number 
of the outer cylinder is Nuo = Nuirj. The effect of the transition between the 
conductive and convective regime on heat transfer is shown in figure|^d). Large 
cr and short F promote convective profiles and thus efficient heat transfer. 




(c) 


(d) 


Figure 3: (color online) (a) Colormaps of the temperature for 77 = 0.5, Re = 30, 
G = 2000, r = 10 and different values of cr. They are plotted in a meridional 
section (r, z) with the inner cylinder on the left hand side. Note that the aspect- 
ratio has been scaled in these plots to enhance the visualization of the contours. 
(b)-(c) Temperature profiles at mid-height and mid-gap, respectively.(d) Vari¬ 
ation of Nu with a for several values of T. The case of cr = 0 has a purely 
conductive temperature profile as in the axially periodic system. 
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3.2. Transition from conductive to convective: the Rayleigh effect 


We quantify the departure of the radial temperature profiles from the con¬ 
ductive regime with the parameter 

fr^iTinid — T’cond)^ 


D = ■ 


i::n 


'2 

cond 


(7) 


where Tmid is the radial temperature profile at mid-height and Tcond denotes 


the conductive profile (2c I of the infinite-cylinder case. 


Figure (a) shows the transition between the conductive and convective 
regimes as a function of the Rayleigh number {Ra = a G) and F. The critical 
Rayleigh numbers Roc have been computed by fixing a (using different values) 
and F, and increasing G until D > 10“^ was reached. For low values of F there 
are significant differences between Rqc calculated with different cr. However, for 
F > 12.5 the transition takes place at nearly the same Rqc regardless of a. This 


interesting behaviour can be explained by analysing equation (2c I. The laminar 


state is steady and axisymmetric, so departures from the conductive regime are 
due to either UbdrTf, or WbdzTi,. The first term is only affected by changes in 
cr and plays a significant role in setups with a relatively short aspect-ratio, as 
the radial velocities arising near the end walls are more vigorous. Consequently, 
variations of Ra^ with a are observed in short setups. 

In contrast, the second term depends on both a and G, but only through 
their product: the Rayleigh number Ra. This becomes progressively dominant 
as F increases and the transition to the conductive regime becomes dependent 
on Ra only. Interestingly, Rac varies linearly with F according to 


Rae(r) = 338.04r - 216.24 


( 8 ) 


This linear behaviour was also observed for fixed ct = 1 by Thomas and De 
Vhal Davis [12], who studied the transition between both temperature regimes 
for natural convection. 

The linear relationship Q separates the regions of parameter space in which 
the temperature regime is conductive and conductive. It provides a criterion 
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(a) 




Figure 4: (color online) (a) Transition from conductive to convective tempera¬ 
ture regime (D > 10“^) as a function of the aspect-ratio F and Rayleigh number 
Ra for r] = 0.5 and Re = 30. (&) Variation of the slope a{rii) in the linear equation 
§ marking the transition between the conductive and the convective regimes. 


to determine a priori whether simulations using axially periodic boundary con¬ 
ditions render a good approximation of laboratory flows. Note however that 


since the temperature profile of the laminar conductive state (2c I depends on 
? 7 , criterion ([^ is valid only for r] = 0.5. In order to investigate the effect of 
curvature in the transition from conductive to convective heat transfer, we per¬ 
formed simulations as those shown in hgure (a) for different values of r/. We 
found that for all investigated curvatures rj G [0.1,0.99] a linear relationship of 
the form 

Rac{r], F) = a(? 7 )r -f 6 ( 77 ) (9) 


is always satisfied, but the slope a depends significantly on rj (see figure]^ (&)). 
The curve 0 ( 77 ) is well approximated by the cubic function 


0 ( 77 ) = -307.277^ -p 750.177^ - 3 I 8.777 + 339.6 (10) 


The offset of the linear equation satisfies b G [—300, —100] for the whole 77 -range 
and so can be neglected in comparison with a(? 7 )r for F > 15. 
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(a) 


(b) 




Figure 5: (color online) (a) Critical stability boundaries {Gc vs Re). The num¬ 
bers that appear on top of the critical curves are the azimuthal modes associated 
to the critical perturbations. (Blue) Solid line and squares are used in the case 
of finite-length cylinders, whereas (green) dashed line and circles stand for the 
marginal stability in the case of infinite long cylinders. The working fluid is air 
(cr = 0.71) and the geometric parameters are rj = 0.5 and F = 10 (in the finite 
case). (5) Detail of the marginal curve for low G in the inhnite cylinders case. 


4. Primary instability 


4.1. Axially periodic disturbances 

The dashed line in figure (a) shows the stability curve of the conductive 
laminar flow to axially periodic disturbances for 77 = 0.5 and a = 0.71. Three 
regions with qualitatively different behaviour can be distinguished. For slow 
or no rotation the laminar flow is destabilised as Ra increases. Buoyancy is 
the indirect force driving the instability, which occurs because the axial veloc¬ 


ity profile (2b) has an inflection point. At low Ra viscosity can stabilise the 


laminar flow, but as Ra increases the magnitude of the axial velocity grows 
proportionally to Ra until the flow becomes unstable at Rac = 5120. The in¬ 
stability pattern consists of axisymmetric convection rolls that travel axially in 
the upward direction [21]. This axial drift of the rolls is caused by curvature. 
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which makes the axial velocity profile of the conductive basic flow radially asym¬ 
metric (see Fig. [^. In the limit of vanishing curvature 77 —> 1 (vertical plates) 
the roll-pattern is steady because the system is left-right symmetric about the 
mid-gap ^5] , 

As the inner-cylinder rotation speed increases the centrifugal force starts to 
compete with buoyancy and the laminar flow becomes first unstable to nonax- 
isymmetric spiral modes. Although centrifugal instabilities are typically axisym- 
metric (i.e Taylor vortices), it is well known that in the presence of a sufficiently 
large axial velocity nonaxisymmetric spiral modes are dominant. This was first 
shown by Snyder [2S] in Taylor-Couette flow with a superimposed axial pressure 
gradient. In our system, as Re increases the critical Rayleigh number for the 
onset of these spiral modes and their wavenumber (n) progressively decrease. 
At {Re,Ra) = (69.42,194.53) the critical curve for the spiral modes intersects 
the curve corresponding to the pure centrifugal instability (i.e Taylor vortices 
with n = 0). Figure ( 6 ) shows a close up to the vicinity of this point. 

4-2. Axially bounded cylinders 

Following Kuo and Ball [15] we chose F = 10. The computed stability 
curve is shown as a solid line in figure (a) and also features an intersec¬ 
tion point around which spiral modes and axisymmetric Taylor vortices com¬ 
pete. In comparison to the periodic case, this intersection point is displaced 
to {Re,Ra) = (66,257.02) because of the hnite aspect-ratio. Here, as in the 
periodic case, and long aspect-ratio experiments of Snyder and Karlson |5|, who 
used F = 337, there is a region in which spiral modes appear at lower values of 
Re than Taylor vortices. However, this region is limited by a turning point and 
when Ra is further increased instabilities are suppressed. The Ekman layers 
cause a slow transition to the convective basic state and enhance heat trans¬ 
fer without the need for instabilities. Overall the stability curves for the finite 
and periodic system are in qualitative agreement up to Ra « 3000. This is in 
line with criterion Q, which predicts the transition between conductive and 
convective states at Ra = 3164. 
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5. Heat transfer for increasing rotation speed 


The stability analysis of the previous section cannot predict the dynamics 
and heat transfer properties of the flow. For these purposes we performed direct 
numerical simulations of the Navier-Stokes equations. We first investigated 
Ra = 1420, so that the temperature for F = 10 is conductive. In the periodic 
case we choose a cylinder length of = 7.48, which allows some freedom in 
the selection of the wavenumber of the vortex-pairs. 

Figure shows the sequence of emerging patterns in the periodic (upper 
row) and finite (lower row) cases as the Reynolds number is increased. For the 
periodic system the primary instability occurs at Re « 39, leading to spiral 
flow with azimuthal wavenumber n = 3 (figure (o))- By further increasing 
Re the axial flow generated by buoyancy loses importance in comparison to 
the centrifugal force. This results in secondary transitions towards spiral flows 
with decreasing azimuthal mode n. The transition from n = 3 to n = 2 takes 
place at Re 52, whereas the transition to n = 1 occurs at Re « 135. The flow 
patterns with n = 2 and n = 1 are shown in figure]^ (6) and (c), respectively. At 
Re « 260 the spiral flow is modulated by a low frequency and the flow becomes 
quasi-periodic. The resulting state preserves the spiral structure with n = 1 
and is characterised by the appearance of spatio-temporal defects (figurej^ (d)). 
Finally, for Re > 280 the centrifugal force completely dominates and spiral flow 
is superseded by Taylor vortices (figure]^ (e)). These results are consistent with 
the physical mechanisms discussed in the linear stability analysis of §4.1| 

The same sequence of flow states was found in the finite case, but with 
slightly different Reynolds number for the transitions. The basic flow loses 
stability to spiral flow with n = 3 for i?e « 50 (figure (/)), whereas the 
transitions to n = 2 and n = 1 occur at Re « 63 and Re « 90 respectively 
(figure {g)-{h))). The transition between spiral flow and Taylor vortices is 
also characterised by the presence of defects, which are much more pronounced 
than in the infinite cylinders case (figure]^ (i)). The onset of Taylor vortices 
(figure]^ (j)) takes place at Re « 130. 
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Figure 6: Evolution of an isosurface of the axial velocity in the periodic (upper 
row) and finite (lower row) systems as Re is varied and Ra = 1420, cr = 0.71, 
T] = 0.5,r = 10 (finite case) and Lz = 7.48 (infinite case). 
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Figure 7: Variation of Nu with Re for Ra = 1420, cr = 0.71 and rj = 0.5 in the 
finite (F = 10) and infinite {L^ = 7.48) cases. 

In figurej^we show the heat transfer rate, normalised by the conductive rate 
as in equation ([^, as a function of increasing Re. As the basic flow is conductive 
for periodic and finite boundary conditions, Nu is in both cases independent of 
Re before the onset of instabilities {Re < 50). In the infinite case heat transfer 
is purely conductive {Nu = 1), whereas in the finite case, Nu is slightly higher 
{Nu = 1.07) due to the heat transfer at the Ekman boundary layers. The onset 
of instability enhances the convective heat transfer significantly and the Nusselt 
number follows a power-law scaling Nu = ARe^, with {A, B) = (0.14, 0.51) for 
the finite case and {A,B) = (0.13,0.53) for the periodic case. 

6. Heat transfer for increasing temperature difference 

We here investigate the dynamics and heat transfer for Re = 50 and in¬ 
creasing Ra. The upper row of figure illustrates the bifurcation scenario for 
axially periodic boundary conditions. Colormaps of the temperature in longi¬ 
tudinal sections at r = are depicted. The onset of instability occurs at 

Ra « 942.88 resulting in a spiral flow pattern with azimuthal mode number 
n = 2. This is stable only in the vicinity of the critical point. A small increase 
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in AT leads to a new spiral state with n = 3 (figure]^ (a)), which remains stable 
for 1065 ;< Ra < 2485. Complex spatio-temporal dynamics emerges as i?a is 
further increased (Ra > 2485). The sequence is as follows. First, spiral states 
turn into wavy spiral flow patterns shown in figures [^(5) and (c), corresponding 
to Ra = 2840 and Ra = 3550 respectively. Note that the dominant spiral mode 
changes back from n = 3 to n = 2 as Ra is increased. These flow patterns are 
similar to those reported in both experimental [S] and numerical studies |27j for 
large F and 77 « 0.8. Subsequently increasing Ra above the wavy spiral flow 
regime the flow becomes chaotic (figure]^ (d)). 

The sequence of flow states obtained in the finite case (lower row of figure]^ 
is qualitatively similar to the periodic case up to the transition to the convective 
basic state, which is found to occur at Ra « 4000, in qualitative agreement 
with criterion (|^. Nevertheless, there are quantitative differences. Here the 
conductive basic flow (figure (e)) remains stable up to Ra « 1420, where 
spiral flow with n = 3 develops (figure (/)). As the radial heating increases 
the dominant spiral mode changes to n = 4 at Ra « 3550 (figure]^ (g)). This 
increase in n with Ra is driven by the stronger axial flow acting on the fluid 
because of thermal buoyancy m- Further increasing Ra leads to the thinning of 
the thermal boundary layers at the sidewalls and the flow transitions gradually 
to the convective laminar flow (figure]^ W). A similar transition from spiral 
flow to the convective state as AT is increased was experimentally observed 
by [13] using an apparatus with F = 31.5. We refer the reader to Ref. [T3] for a 
detailed description of these instabilities. 

The evolution of the normalised Nu with Ra clearly illustrates the different 
dynamical behaviours in the finite (figure]^ (a)) and infinite (figure]^ ( 6 )) sys¬ 
tems. The largest differences occur for weak rotation {Re = 20). In the infinite 
case , Nu = 1 remains constant up to Ra « 7810, indicating that heat trans¬ 
port is exclusively conductive up to reach the onset of instability. In contrast, 
in the finite case, the Ekman layers provide an efficient mechanism to transfer 
heat, which results in an almost linear increase of Nu with Re. Such growth 
is also observed for spiral flow. Re = 50 and Re = 70. When Ra exceeds the 


17 



Figure 8: (color online) Evolution of the temperature in the periodic (upper row) 
and finite (lower row) systems as Ra changes and Re = 50, cr = 0.71, rj = 0.5, 
r = 10 (finite case) and = 7.48 (infinite case). Colormaps are plotted in 
longitudinal sections (0, z) at r = , with the exception of the axisymmetric 

states (e) and (h), which are plotted in a meridional section (r, z). Note that 
these plots have been resized and do not correspond to the aspect-ratio of the 
investigated system. 
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Figure 9: Variation of Nu with Ra as a function of Re for a = 0.71 and rj = 0.5 
in the (a) finite (F = 10) and (6) infinite = 7.48) systems. 

value corresponding to the upper part of the marginal curve in figure (a) the 
flow becomes again axisymmetric and Nu solely depends on Ra. As a result, 
curves corresponding to different values of Re collapse. In the infinite case, 
the behaviour of Nu within the region of spiral flow. Re = 50 and Re = 70, 
reflects the transitions described above. The spiral flow {Ra < 2750) resulting 
from the primary transition provides efficient heat transfer, leading to a rapid 
linear growth of Nu as Ra is increased. In this flow regime the values of 
in the infinite case are slightly larger than those in the finite case. The tran¬ 
sition towards wavy spiral flow results in a sudden decrease of the convective 
heat transport, consistently with the results of HU. After this initial drop, fur¬ 
ther increasing Ra towards fully developed turbulent flow is accompanied by a 
progressive growth in Nu., which is however significantly lower than the linear 
growth from the spiral flow regime. 

7 . Discussion 

Heat transfer in laminar flows between two concentric vertical cylinders can 
be of conductive or convective nature. Regardless of the curvature of the cylin- 
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ders, parametrised here by the radius ratio 77 , the transition between these 
regimes is governed by the Rayleigh number Ra and the length-to-gap aspect 
ratio r. In particular, for a given T there is a critical Ra above which the lam¬ 
inar state is convective. We here determined numerically a linear criterion ([^, 
which applies to all geometries and is independent of the Prandtl number as 
long as r > 15. For smaller F the transition to the convective state is influenced 
by a because of the stronger effect of the Ekman layers. Although the criterion 
was obtained for a fixed Reynolds number Re = 30, we provided evidence that 
it appears valid as Re increases. This was done by comparing results obtained 
with axially periodic boundary conditions, for which the laminar flow is always 
conductive, to results obtained with physical (no-slip) boundary conditions. 

If the laminar flow is conductive and Re is increased while Ra is kept fixed, 
our direct numerical simulations show that the transition to chaotic flow is very 
similar in the periodic and finite cases. This was shown by using 77 = 0.5, 
a — 0.71, Ra = 1420 and F = 10 in the finite case, following Kuo and Ball [T5] . 
All their reported transitions were reproduced here, however, we found a notable 
difference in the onset of spiral flow, which occurs for lower Re in our simulations 
and corresponds to a critical spiral mode n = 3 instead of tt, = 2. This can be 
explained because Kuo and Ball neglected the correction of the velocity field 
in their simulations, so their velocity field was not mass conserving. Omitting 
the correction in our simulations we reproduced their results and found that its 
relevance diminishes as Re increases, so that the secondary transitions reported 
here coincide with those in m- 

We further repeated these simulations by using periodic boundary conditions 
in a domain of length = 7.48. We found that the onset of instability occurs 
at slightly lower Re than in the finite case, where the Ekman layers act as 
frictional layers, slowing down the internal flow and thus delaying the primary 
and secondary instabilities. In contrast, as the Re is further increased the 
remaining instabilities, culminating in Taylor-vortex flow, are favoured by the 
presence of axial end walls. In this case the Ekman circulations at the end wall 
enhance the Taylor vortices, as happens in the isothermal case. 
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The heat transfer, quantified by the convective Nusselt number ([^, follows 
a simple power law {Nu = and Nu = in the finite and 

infinite cases respectively). The exponents are in good agreement with the 
analytical correlation provided by Dorfman [^. A similar exponent {B = 0.47) 
has also been found by Viazzo and Poncet m who numerically explored a 
region of parameter space similar to that in figure but using an apparatus 
with 77 = 0.8 and T = 80. The small discrepancies between the correlations 
obtained in the finite and infinite cases occurs because of the Ekman layers, 
whose importance grows as T is reduced. 

If Re is kept fixed while Ra increases, the transition to chaotic flow and 
heat transfer of the periodic and the finite system remain similar as long as 
criterion is approximately satisfied. In the finite system, further increases in 
Ra result in a transition to the laminar convective state, whereas the periodic 
system becomes gradually more turbulent. Thus there are two counter-acting 
mechanisms as Ra increases. Instabilities and turbulence are stimulated be¬ 
cause of thermal driving. However, thermal driving results in a thinning of the 
thermal boundary layers because of end-wall effects. This quenches turbulence 
and causes the transition to the laminar convective state. Noteworthy, at a fixed 
Ra the laminar convective state is much more efficient in transferring heat than 
the turbulent flow of the periodic system. 

8. Conclusions 

The role of axial boundaries in laterally heated and rotating flows was in¬ 
vestigated with linear stability analyses and direct numerical simulations. Both 
periodic and no-slip boundary conditions were considered in the axial direction. 
The main results can be summarised as follows: 

• Criterion ([^ determines whether laminar heat transfer in an axially bounded 
system is of conductive or convective nature. It is independent of Prandtl 
number and robust with Reynolds number as long as P > 15. 
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• If the laminar state is conductive axially periodic boundary conditions 
correctly describe the dynamics and heat transfer in laboratory setups. If 
the laminar state is convective the infinite-cylinder approximation may be 
used but the basic flow must be modified to include axial stratification [53] . 

• A variety of flow patterns was observed. The underlying flow instabil¬ 
ities are driven by the interaction of two physical mechanisms: thermal 
buoyancy generates an inflectional axial velocity profile, whereas the inner- 
cylinder rotation promotes centrifugal instability. 
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